Self-organized criticality induced by quenched disorder: experiments on flux 

avalanches in NbH a . films. 
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We present an experimental study of the influence of quenched disorder on the distribution of flux 
avalanches in type-II superconductors. In the presence of much quenched disorder, the avalanche 
sizes are power-law distributed and show finite size scaling, as expected from self-organized criticality 
(SOC) . Furthermore, the shape of the avalanches is observed to be fractal. In the absence of quenched 
disorder, a preferred size of avalanches is observed and avalanches are smooth. These observations 
indicate that a certain minimum amount of disorder is necessary for SOC behavior. We relate these 
findings to the appearance or non-appearance of SOC in other experimental systems, particularly 
piles of sand. 
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Self-organized criticality (SOC) [1] has generated great 
interest over the last 15 years mainly due to its wide 
range of possible applications in many non-equilibrium 
systems [2,3]. However, progress has been hampered by 
the fact that clear, tell-tale signs of criticality, such as 
finite-size scaling in the distribution of avalanches, have 
only been observed in very few controlled experiments 
[4]. Recently however, there have been a number of ex- 
perimental observations of criticality in both two- [5,6] 
and three dimensional systems [7-9]. However, the criti- 
cal ingredients to obtain SOC in an experimental system 
still remain obscure. 

Recent theoretical advancements have studied the na- 
ture of the criticality in SOC and made a link with phase 
transitions describing how a moving object comes to rest 
[10,11]. Such absorbing state phase transitions are closely 
related to the roughening of an elastic membrane in a 
random medium [12]. In these theoretical models, there 
needs to be an absorbing state phase transition underly- 
ing the process at hand in order to observe SOC. This 
critical state is reached by a self-organization process 
[13], which depends on slowly driving the system away 
from its state where everything is at rest. If the driv- 
ing is not slow enough, the relaxations to the critical 
point may be disturbed by the driving, such that no crit- 
ical state is reached [14]. Too strong driving may have 
occurred in some of the early experiments, particularly 
those performed in rotating drums and where the grains 
were dropped from a considerable height. The absence of 
an underlying phase transition, however, would be more 
fundamental and detrimental to SOC. In the case of ab- 
sorbing state phase transitions, it is known that the pres- 
ence of quenched disorder plays an important role in the 
nature of the critical point, such that this may also be 
an important ingredient in obtaining SOC behavior. 

Here we present an experimental investigation of the in- 
fluence of disorder on the appearance or non-appearance 
of SOC. Therefore it is necessary to have a system, where 
the quenched disorder can be experimentally changed 
while leaving all other aspects the same, as well as having 



a system which has been shown to show SOC at least in 
some circumstances. We study the avalanche dynamics of 
magnetic vortices in the type-II superconductor Niobium 
in the presence of Hydrogen impurities using magneto- 
optics. As first noted by de Gennes [15], the penetration 
of a slowly ramped magnetic field into a type-II supercon- 
ductor has a strong analogy to the growing of a sand pile, 
the archetypal example of SOC. It has been shown in the 
past, that these vortex avalanches are distributed accord- 
ing to a power-law [16] and more recently that they obey 
finite-size scaling [9]. However, due to the presence of 
pinning in the system, studying magnetic vortices allows 
a detailed investigation of the influence of quenched dis- 
order on their avalanche distribution and structure. Nb 
is particularly suited for this purpose, as it easily takes 
up H impurities (which locally destroy superconductiv- 
ity, see below) [17], which can be added to the sample via 
a contact gas, thus allowing an experimental control of 
the amount of quenched disorder in the system. Further- 
more, studying magnetic vortices has the advantage that 
kinetic effects in their dynamics, which are thought to 
have hampered some sand-pile experiments [5], are nat- 
urally absent due to the over-damped dynamics of the 
vortices [18]. 

The experiments described here were carried out on 
a Nb film of a thickness of 500 nm, evaporated under 
ultra high vacuum conditions on an 'R-plane' sapphire 
substrate. The films were then covered with a 10 nm Pd 
caplayer in order to facilitate the catalytic uptake of H 
into the film [19]. 

The local magnetic flux density, B Zl just above the 
sample is measured using an advanced magneto-optical 
(MO) setup, directly yielding the local Faraday-angle in 
an Yttrium-Iron Garnet indicator, using a lock-in tech- 
nique [20] . Images are taken with a charge- coupled device 
camera (782 x 582 pixels) with a resolution of 3.4 /im per 
pixel. The sample is placed in a cryo-magnet and cooled 
to 4.2 K in zero applied field. The external field is then 
ramped from to 20 mT in steps of 50 /zT, where the 
flux in the sample is relaxed for 3 seconds before an MO 
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picture is acquired. This sequence is repeated twice for 
each H concentration to check for reproducibility. 

The quenched disorder in the sample is increased af- 
ter such a sequence of two field sweeps by equilibrating 
the sample with a certain higher partial pressure of H 
at room temperature for one hour. The partial loading 
pressures used in the experiments discussed here were 
80, 260, 1130, and 1810 Pa. We estimate that the H im- 
purities present in the as grown sample correspond to a 
partial loading pressure lower than roughly 10 Pa. Af- 
ter equilibration, ensuring a uniform distribution of H 
in the Nb film, the sample is cooled down again. In 
this cooling procedure, phase separation occurs in H-rich 
and H-poor regions [21], where superconductivity is sup- 
pressed in the H-rich phase [22]. Thus the clusters of 
H-rich phase, which are in the order of 0.1 to 1 /zm in di- 
ameter [21], act as effective pinning sites for the vortices 
and thereby introduce quenched disorder. 
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FIG. 1. a) The magnetic flux landscape in the center of 
the Nb film without added Hydrogen impurities. The mag- 
netic flux density in the sample is indicated by the brightness, 
where white corresponds to 32 mT and black corresponds to 
mT (see inset), b) The flux landscape in the same sam- 
ple after the application of a H pressure of 1810 Pa at room 
temperature. As can be seen, the flux structures are much 
more disordered and branched, c) Magnetic flux avalanches 
in the absence of H treatment, as given by the subtraction 
of two consecutive MO images. Note that the magnification 
and the grey-scale is different from a), white corresponds to 
mT and black corresponds to 3 mT. The avalanches have a 
predominant size and are rather smooth and plume-shaped, 
d) Flux avalanches after H treatment with 1810 Pa gas pres- 
sure at room temperature, again obtained from a subtraction 
of consecutive MO images (grey scale and magnification are 
as in c). The avalanches are more fractal and branched than 
in c) and do no longer have a preferred size (see text). 

In this manner, we obtain a collection of magnetic flux 
landscapes as shown in Fig. la), which are then analyzed 
in terms of the avalanches that have occurred between 
time steps. In order to obtain the amount of magnetic 
flux that has been transported in one avalanche, we first 
average each image over a set of 4 x 4 pixels to reduce 
noise. Then two consecutive images are subtracted lead- 



ing to an image like that shown in Fig. lc). As can be 
seen in the figure, in each image we observe a number 
of avalanches. These are identified individually, using a 
threshold of 0.3 mT in AB Z . The amount of flux in each 
of these avalanches is subsequently determined from in- 
tegrating the difference in flux density over the area of 
the avalanche: 



AB z dxdy. 



(1) 



From the spatial resolution of the setup and the thresh- 
old level of the avalanche identification, we can determine 
the smallest avalanches that are still resolved to contain 
about 20 $o, where $o = h/2e is the magnetic flux car- 
ried by each vortex. In order to check for finite size scal- 
ing, we discard avalanches exceeding a linear extent of 
200, 100, and 50 fim from the analysis. This corresponds 
to decreasing the system size accordingly, as the linear 
extent of an avalanche cannot exceed the system size. 

As can be seen in Fig. lc), when the sample has not 
been treated with H, the avalanches have a characteris- 
tic structure and area. This leads to a peak in the size 
distribution, which disappears as avalanches exceeding 
the preferred linear extent are discarded. Thus in the 
absence of H-induced quenched disorder, no SOC behav- 
ior is observed, as there is no finite-size scaling of the 
avalanche distribution. This is shown in Fig. 2a), where 
we show the scaled avalanche size distribution for the 
sample without H treatment. Here, the size distribution 
has been logarithmically binned, such that avalanche size 
steps in the histogram are separated by a constant fac- 
tor rather than a constant step width. Moreover, the 
histograms are scaled with s T , which leads to horizontal 
lines in case the avalanches are power-law distributed. 
The displaced flux per avalanche is scaled by L~ D , where 
L is the scale above which avalanches are discarded and 
D is the fractal dimension of the avalanches. In order to 
obtain the 'collapse' in Fig. 2a), we have used r = 1.35 
and D = 2.75. The value of r is determined from a best 
collapse of the data onto a horizontal line for small s. The 
value of D is found from the best collapse of the cut-off 
values. The nice collapse of cut-off values indicates that 
the avalanches in this case do have a definite fractal di- 
mension, which is close to 3, as one would expect given 
their smooth appearance in Fig. lc). The peak that can 
be observed in the distribution is due to the predomi- 
nance of avalanches containing about ~ 2500 $o ; which 
can also be inferred from Fig. lc). 

Turning to the data in Fig. Id), it can be seen that in 
the presence of strong quenched disorder (i.e. after treat- 
ing the sample with 1810 Pa of H gas), the avalanche 
structure is markedly different from that in the pristine 
case. Now, the flux moves in a much more branched 
way and furthermore a preferred size of avalanches is 
absent. Again, we have checked this with a finite-size 
scaling analysis of the avalanche size distributions, see 
Fig. 2b). The size distribution is again logarithmically 
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binned and scaled in order to obtain a curve collapse. 
Here, in contrast to Fig. 2a), there is good curve col- 
lapse and thus we can determine the size distribution 
exponent and avalanche dimension to be r = 1.07(2) and 
D = 2.25(5), where we observe a power-law distribution 
over more than two orders of magnitude. 
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FIG. 2. a) The flux avalanche size distributions in the ab- 
sence of Hydrogen treatment. Data are shown for sets of 
avalanches not exceeding different linear sizes of L = 200, 100, 
and 50 /im. The different distributions are vertically scaled 
with s T and the sizes are horizontally scaled with L~ D , in 
order to obtain the usual curve collapse for finite size scaling. 
However, as avalanches of ~2500 <I?o are preferred, no good 
curve collapse is obtained in this region. For small avalanches 
and avalanches of the cut-off size, the curves do collapse, such 
that the fractal dimension of the avalanches can be deter- 
mined, b) The avalanche size distributions scaled in the same 
way, for experiments after H treatment at 1810 Pa of the sam- 
ple. Here a good curve collapse is observed, with power-law 
scaling over more than two decades. This is clear evidence 
that in the presence of quenched disorder, SOC is present in 
the vortex avalanches. 

This overall behavior of flux avalanches in the presence 
of changeable quenched disorder has also been studied in 
molecular dynamics simulations of a collection of vor- 
tices [23]. This is a great advantage of the vortex system 
over the study of sand piles, as the interactions between 



vortices are well known [18] and can thus be simulated 
in order to separate collective effects from microscopic 
dynamics. In their study, Olson et al. [23] have found 
that on increasing the pinning density, the avalanche dis- 
tributions change from showing a preferred size (with a 
power-law distribution at small sizes) to a power-law dis- 
tribution. Furthermore, in the case of strong pinning, 
the exponent r obtained in these power-law distributions 
is between 1 and 1.4 in very good agreement with our 
experimental values. Moreover, in molecular dynamics 
simulations, the movements of all vortices can be followed 
and it turns out that in the case of low pinning density, 
the avalanches form channels [23], which is highly remi- 
niscent of the structures seen in Fig. lc). 

Let us now turn to the general implications on the ap- 
pearance of SOC that can be drawn from the above ob- 
servations. As we have shown, in the absence of quenched 
disorder, even the non-kinetic vortices [18] do not show 
finite-size scaling in their avalanche distribution. Fur- 
thermore, it has been shown previously that steel balls in 
a two dimensional pile can show SOC [6] when quenched 
noise is introduced via a fixed random arrangement of 
balls at the base. These two observations together seem 
to indicate that kinetic effects may be less important 
than previously thought [5], while it is the presence of 
quenched noise that leads to the appearance of SOC. 
This is also corroborated by molecular dynamics simu- 
lations on vortices, using over-damped dynamics finding 
the same result [23]. 

The importance of quenched disorder can be best seen 
in the structure of the avalanches, which is intimately 
connected to the surface of a pile or flux landscape [7,24]. 
As the number of pins is increased, the vortex avalanches 
have to accommodate to their presence. This frustra- 
tion leads to a much more branched and open avalanche 
structure. This is shown in Fig. 3, where the circles show 
the values of D for different H loading pressures. As 
can be seen, at low pinning density, the avalanches have 
a dimension close to 3, indicating a smooth structure, 
whereas with increasing H content, the avalanche dimen- 
sion decreases. Due to the more complex structure of an 
avalanche, a bigger range of possible avalanche sizes can 
be explored over the area of the sample. Thus there is 
no longer a preferred size of avalanches and hence SOC 
behavior occurs. Again, we can compare our results to 
those of Altshuler et al. [6], where in the absence of SOC, 
only a small part of the pile took part in avalanches and 
the pile surface overall was much smoother than in the 
case where SOC was observed. 

In order to quantify the roughness [25] of the 
flux landscape, we determine the width, ((B z (x,y) — 
(B z (x,y)}y) 2 y/ 2 : of the flux landscape after subtraction 
of the average profile, (B z (x,y)) y . As can be seen in 
Fig. 3, with increasing H content, the flux landscape be- 
comes more rough as given by a fourfold increase in sur- 
face width. Thus as has been observed in experiments 
using steel balls [6] as well as in rice piles [7,26], SOC be- 
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havior in the flux landscape is coupled to the observation 
of a rough pile surface. 
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FIG. 3. The fractal dimension of the avalanches (cir- 
cles) and the surface width (see text) of the flux landscape 
(squares) as a function of Hydrogen partial pressure. Data 
for the pristine films are shown at 10 Pa, the estimated max- 
imum contamination H pressure. The fractal dimension de- 
creases with increasing quenched disorder, corresponding to a 
more branched and rough structure, as seen in Fig. Id). The 
surface width increases with quenched disorder, which indi- 
cates that the flux landscape becomes increasingly rough as 
H is introduced into the sample, as can be seen in Fig. lb). 

In conclusion, we have shown in the case of high pin- 
ning density that the magnetic flux avalanches observed 
in a Nb film show finite-size scaling, which implies the 
presence of SOC in the system. However, in the presence 
of a low pinning density, a preferred size of avalanches is 
found and no finite-size scaling is observed. This demon- 
strates the importance of quenched disorder in the system 
in order to obtain SOC behavior. The exact position of 
the transition from non-SOC to SOC behavior cannot be 
determined very accurately due to experimental limita- 
tions (e.g. 'noise' in Fig. 2 and limited accessible range of 
L), but might also be intrinsically smooth if one considers 
the gradual changes observed in Fig. 3 for the avalanche 
dimension D and surface width. Our experimental ob- 
servation of this transition as a function of disorder is 
in agreement with previous numerical work on cellular 
automata [27] and with molecular dynamics simulations 
of magnetic vortices [23]. In addition, our findings are 
consistent with the point of view of SOC as an absorbing 
state phase transition [10,13]. In absorbing state phase 
transitions, such as directed percolation [25], the density 
of quenched disorder is a critical parameter, whereby a 
phase transition can be induced. The presence of such 
an underlying phase transition is a necessary condition 
in order to obtain SOC in the model of Vespigniani et al. 
[13]. Thus again, this view advocates that the increase 
of quenched noise can lead to the appearance of SOC, as 
indeed we have found experimentally. 
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